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ABSTRACT 

We examine nonlinear oscillations of slender tori in the vicinity of black holes and compact stars. These tori represent useful probes 
of the complicated, nonlinear dynamics of real accretion disks and provide at least qualitative understanding of their oscillations. We 
demonstrate that epicyclic modes of such tori are weakly coupled due to the pressure and gravitational forces. We explore all possible 
resonances between two epicyclic modes up to the fourth order. We show that the strongest resonance between axisymmetric modes is 
3:2. In addition, any resonance between an axisymmetric and a non-axisymmetric mode is excluded due to axial and equatorial-plane 
symmetries of the equilibrium torus. We examine a parametric excitation of vertical axisymmetric oscillations by radial oscillations 
in the 3:2 resonance. We show that the resonance may be significant only for high-amplitude radial oscillations. 

Key words, accretion, accretion-discs - oscillations - black hole physics - quasi-periodic oscillations 



1. Introduction 

In the past thirty years, the oscillations of fluid tori orbiting 
around massive compact objects were studied systematically, 
particularly in the context of the stability of thick accretion disks 
dPapaloizou & Prinrielll984lll985UBlaeslll985l;[Goldreich et al.l 
1986; Goodman et 19871 iNarayan et al.ll 19871 and many otiV 
ers). Torus oscillations have attracted renewed attention fol- 
lowing the discovery of kilohertz double-peak quasi-periodic 
oscillations (QPOs ) in the light curves of sev eral accreting 
black holes (see iMcClintock & Remillardl 120061 for a review). 
Frequencies of these oscillations are often in the 3:2 ratio and 
they are detected mainly when a source is in a so-called steep 
power-law state. The geometry of the accretion flow in this state 
is not yet clear, nevertheless a fluid torus created by substan- 
tial pressure gradients in the inner parts of the accretion disk 
and supporting discrete trapped oscillations is a plausible pos- 
sibility. Such tori appear in many non-radiative global simula- 
tions as natural components of magnetohydrodynamic turbulent 
accretion flows dBalbus & Hawlevl2002llDe Villiers et alj2003t 



iMachida et aT] |2006) 

In some models, QPOs are identifie d with two d i fferen t 
modes of torus oscillations. For example. iRezzolla et al.l d2003l) 
described QPOs as the two lowest-order p-mode oscillations 
of a polytropic torus of small thic kness with constan t angular- 
momentum distribution. Similarly, iBlaes et ail (|20p6b explored 
linear oscillations in slender tori and identified the twin-peak 
QPOs with the vertical epicyclic mode and the breathing mode. 
Although frequencies of these modes depend on the position of 
the torus in the accretion flow, their ratio remains close to 3:2 in 
a wide range of radii. 

In this context, |Xb ramo wicz & Kluzniak! (|2001) proposed 
an interesting general idea that a nonlinear resonance between 
two modes of accretion-disk oscillations is responsible for ob- 
served QPOs in both black-hole and neutron-star sources. Being 
more specific, they identified them with the radial and verti- 
cal epicyclic modes of accretion tori dKluzniak & Abramowiczl 
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2002) and, adopting the approximation of a slender torus, they 
showed that the existence of t hese modes is quite inde pendent 
of the torus equation of state (Abram owicz et al1l2006l) . Later, 
Blae s et aT] d2007l) demonstrated that the epicyclic modes also 
exist in thicker polytropic tori and found approximate expres- 
sions for their eigenfunctions and eigenfrequencies. 

The resonance arises due to nonlinear coupling be- 
tween the epicyclic modes. I n the scenario suggested by 
Kluz niak & Abramowiczl d2002l) . one mode is first excited ex- 
ternally, either by a n external periodic forcing due to rotation of 
the central object dLee et alj|2004t lLeeH2005l) or by a stochasti c 
forcing due to turbulence in the accretion flow (IVio et al.l l2006). 
The other mode is then continously supplied by the energy from 
the former one via the parametric resonance. The nonlinear in- 
teraction may be strong enough and the amplitude of the second 
mode may eventually exceed the amplitude of the first mode. In 
the strong gravitational field, the strongest resonance occur when 
the frequency ratio of the vertical and radial oscillations is close 
to 3:2. 

In the context of QPOs, the reso nance effects in thin ac- 
cretion disk were already studied by iKatol d2003l 12004 120081. 
and references therein). In these models, QPOs are identified 
with waves (corresponding to either g-mode or p-mode oscil- 
lations), parametrically excited by the deformation of the disk 
(warp or precession). Mutual nonlinear interactions between 
different modes of a th in disk were recently studied also by 
Fogelstro m et all d2008) using a local approximation. 

In this note, we further examine the importance of nonlin- 
ear coupling between epicyclic oscillations of the slender torus. 
Since the geodesic equations are separable, earlier studies, based 
on epicyclic motion of test pa rticles, used an additional 'ad-hoc' 
force to initiate the resonance (Abramo wicz et al.l2003 ; Rebusco 
2004). We show that in the case of epicyclic modes of fluid 
torus, the situation differs and both epicyclic oscillations are nat- 
urally coupled due to pressure and gravity. We also estimate the 
strength of this coupling. 

The plan of the paper is following. In Sect. [2] we briefly in- 
troduce the approximation of slender tori and review the for- 
malism that we use to calculate nonlinear interactions among 
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slender torus modes. This formalism is similar to that used fre- 
quently for modeling nonlinear oscillatio ns of rotating stars (e.g. 
Dyson & Schutzl 119791: ISchutzl Il980allbi: iKumar & Gold reich 
1989: IWu & Goldreichl 120011: ISchenk et alJ l2002t lArras et alJ 
2003). Section[3]describes the coupling coefficients. We estimate 
the relative importance of pressure and gravitational coupling in 
different types of torus oscillations. In Sect. [4] we discuss pos- 
sible resonances between epicyclic modes based on their sym- 
metry properties, and describe a particular example of the 3:2 
epicyclic resonance. Sections [5] and [6] are devoted to discussion 
and our conclusions. 



2. Formalism 

A stationary configuration of the Newtonian slender torus in a 
general axi s ymme tric gravitational field <I>(r, z) is described by 
Blae s et all d2007l) . The torus consists of a polytropic fluid with 
the equation of state p oc p 1+1 ^", where p and p are the pressure 
and the mass-density at a given point, and n is the polytropic 
index. In cylindrical coordinates \r, <p,z), the velocity of the sta- 
tionary unperturbed flow is purely azimuthal v = Qc^ (e^ is the 
unit vector in the azimuthal direction) and constant at cylinders 
r = const, i.e. O = O(r). Furthermore, the torus is symmetric 
with respect to the equatorial plane. The density and pressure 
profiles of the torus can be expressed as p(x) = pof"(x) and 
p(x) = pof n+l (x), where po and po are the density and pressure 
at the center of the torus (the circle r — ro; hereafter the subscript 
'0' refers to an evaluation at this point) and f(x) is an auxiliary 
function (/ = and 1 correspond to the boundary and the cen- 
ter of the torus). The size of the torus depends on the slender- 
ness parameter yS, which is defined by /3 2 = 2(n + l)po/(po*Q&o)' 
The azimuthal velocity at the torus center is given by the local 
Keplerian angular frequency Qo = (r 3 r <D)o. 

In the limit of (3 — > 0, the torus with a constant-angular- 
momentum distribution has an elliptical cross-section that have 
lengths of the major and minor axes in the ratio of the radial and 
vertical epicyclic frequencies measured at the center of the torus, 
u 2 - (^3>)o+3Qq and u 2 - (d?0)o. The distance from the outer 
or inner edge to the center of the torus is Ar = p>o/o> r + 0(J3 2 ), 
where a) r = cj r /Qo. 



2. 1 . Linear perturbation 

Since the equilibrium configuration is stationary and axisym- 
metric, all linear perturbations are proportional to exp[-i(wf - 
nuf))], where co e C and m e Z are the eigenfrequency and 
the azimuthal wavenumber of the perturbation. The eigenfunc- 
tions of the torus are tradit ionally expressed in terms of the 
Papaloizou & Pringle (1984) variable W = -dp/(crp), where 
6p is the Eulerian pressure perturbation and cr = o> - mO is 
the eigenfrequency measured in the system comoving with the 
equilibrium flow. In next, these functions are referred to as the 
Eulerian eigenfunctions. 

Alternatively, the same perturbations can be expressed in 
terms of the Lagrangian displacement 



r mk 2 W -2cr£lrd r W i(mcrW - 2Q.rd r W) 



2Q.r(cr 2 - k 2 ) 



(cr 2 — k 2 )r 



d z W 
cr 



(1) 



where k 2 - r Q 3 d{ 2 /dr and l(r) is the angular-momentum dis- 
tribution, £(r) = r^Qir). The vector field £ is referred to as the 
Lagrangian eigenfunctiorfl. 

For infinitely slender torus, the Eulerian eigenfunctions of 
the radial and vertical epic yclic modes are W r oc (r- ro) exp(im<p) 
and W v oc zexp(im0) (see Bl aes et aTll2006l) . Each of the eigen- 
functions corresponds to both a positive and negative corotation 
eigenfrequency of the mode, <r = ±u> r and cr = ±0i z , re- 
spectively. In the Lagrangian description, the eigenfunctions and 
eigenfrequencies of the epicyclic modes are given by 



CJ r r 



w r = mOo ± cl>i- 



and 



£ v =yt v {0, 0, l}e 



(2) 



(3) 



where & T and Jl v are normalization real constants. The eigen- 
functions of the two modes with opposite eigenfrequencies are 
complex conjugated, i.e. if (u>, g) represents the solution to the 
eigenvalue problem then another solution is {-u, £*). 

2.2. Nonlinear coupling 

The nonlinear oscillations of polytropic torus are governed by 
the partial differential equation 



D 2 £ ! 



-V; [(y - l)p(V • &gfi + pV l e] + &VjV<b = (4) 



where D/Df - d, + v^V^ is the Lagrangian flow derivative, 
y = 1 + l/n is another polytropic index and g'-i denotes con- 
travariant components of the metric tensor. All nonlinear terms 
are absorbed into the acceleration term a' n (g) on the right-hand 
side of the equation. 

At a given time, the Lagrangian displacement of the nonlin- 
ear oscillations can be expressed as a linear combination of the 
Lagrangian eigenfunctions £J3, 



£(r,*) = 2 c A«fi«, 



(5) 



where the capital Latin indices symbolize different modes. The 
coefficients ca describe the time-dependent instantaneous am- 
plitudes and phases of different linear modes. When this form is 
substituted into Eq. (0J, we obtain a system of many oscillators 



dcA i , 

— — +ioj a ca = —T A {ci). 
at bA 



(6) 



Different oscillators are coupled by nonlinear functions Ta(ci). 
In the absence of nonlinearities, the solutions of Eq. © are har- 
monic functions c A ~ exp(aM0 an d Eq. (O becomes just a su- 
perposition of linear modes. The coefficients b A depend solely 
on £4 and takes into account the normalizations of the modal 
eigen functions. They are given by dSchutzlll980at ISchenk et all 
2002), 



= 2jTp{ l7 - x [(£) 



2 + (S) 2 + (^) 2 ]-2irQ£^}dV, (7) 



1 The differences betwee n equation |Q3 and the equation (2.14) of 
Papaloizou & Pringle | |1985|) are due to a different definition of the 
corotation frequency cr 

2 In the case of Jordan chain modes the expansion |T5} must be com - 
pleted by the associated Jordan-chain vectors (see Schenk et al. 2002). 
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where V is the torus volume. The nonlinear coupling functions 
can be expanded into a Taylor series, 



Ta = ^ «abc cbCc + Y] kabcd c b c c c d + 

B,C 

Z 



and 

ABCDE 



KABCDE CbCcCdCe 



(i + 6y - 4y + 3y ^AmncnDm + 



lOyCy - 3)Ti (A riBT]criDE) + 20y7] (A i]BTlCDE) + 



(8) 



15(y - l)T]( A r]BcriDE) + 2Qj] (A Br]CDE) dV, 



B,C,D,E 



The coefficients of this expansion kabc, kabcd and kabcde are 
the three-mode, four-mode and five-mode coupling coefficients. 
These are studied in detail in the following section. 

Equation (O is not applicable to Jordan-chain modes, for 
which bA = 0. This is not however the case for the epicyclic 
modes. It can be verified using Eqs. ©, 0) and (0 that 



ABCDE ' 



"4! I' 



p^M^V i V,V z V m V n 0)dV > 



(16) 



(17) 



fo r = ±2w,.M5 ? [ 2 , b v = ±2uj z M3l;, 



where 



M 



I 



p(x)dV : 



27r 2 /j 2 ^ P() 
(n + 1) u) r u> z 



(9) 



(10) 



is the torus mass. In the limit of an infinitely slender torus, 
bA — for corotation modes. They correspond to Jordan chains 
of length 1 and become unst able when ft > 0, developing the 
Papaloizou-Prin gle instability (|Blaesll985l) . This agrees with the 
general theory o f Schutz( 19803jT~ 

The decomposition of the nonlinear solutions into eigenfunc- 
tions of linear modes is a common procedure in the theory of 
stellar oscillations. Up to the second order in coefficients ca, our 
Eqs. (H), © and dp are identical to Eqs. (4.12) and (4.13) of 
ISchenk et alj (|2002), who considered only three-mode nonlinear 
interactions. 

In principle, a similar approach can be applied to a more gen- 
eral case of MHD flow. An equation that models the nonlinear 
evolution of the oscillation modes, would have a similar form 
as Eq. d6j; the specific expression for the coupling coefficients 
would, however, differ from the present case. 

3. Coupling coefficients 



where brackets in the indices denote the symmetrization. For a 
simpler notation, we introduce scalars 

t\a = (Vj&), (18) 

rjAB = (V^XVyQ (19) 

TlABC = (V^XV^XV^). (20) 

The three-mode pres s ure-coupling coefficients we re derived 
by iDziembowskH d 19821). iKumar & Goldreichl d!989l) . and re- 
viewed recently by ISchenk et alJ d2002|). The four-mode cou- 
pling coe fficients were deriv ed byjVan Hoolst & Smeyers ( 1993) 
(see also IVan Hoolstl 1 19941) who considered the more general 
case of the self-gravitating isothermal star. 

Before calculating exact values, we review the necessary 
conditions for coupling coefficients to be nonzero and explore 
the importance of different kinds of coupling for different types 
of torus oscillations. For this purpose, it is useful to introduce a 
multi-index notation; the multi-indices are denoted by bold-face 
letters, Xa = Xa^.a,, and their absolute values are given by the 
number of the indices, |A| = n. 

3. 1 . Selection rules 

The necessary conditions for non-zero coupling coefficients fol- 
low from the symmetry properties of the integrands in Eqs. (fT2l> — 
dTTl i. and are natural generalization of the well-known sele ction 
rules for three-mode coupling (see e.g. ISchenk et al.l l2002). The 
azimuthal selection rule states that the integrands cannot depend 
on the azimuthal angle 0, 



Individual fluid elements of the torus move under the combined ^H, 

influence of the gravitational and pressure forces. For this rea- 2_i mAt ^ ^ ^ k a -0. 
son, it is beneficial to separate the contributions of pressure and 
gravity and express coupling coefficients, in general, as 



(21) 



k=\ 



K = K (P) + K (g) . 



(ID 



These two contributions are referred to as the pressure-coupling 
and gravitational-coupling coefficients. Up to the fourth order, 
they are given by 

k abc = J X P ( (r ~ 1)27M7?B77c + 3(r ~ 1 ) J M J 7 2f Q + 

2/7(^c)}dV, (12) 

"aIc = ~\ f v P&U k c WjV^dV, (13) 

k abcd = ~ j[ J v yp {( 3 ~ 3 y + y 1 } iatibicid + 8 noniscD) + 

6(7 - 2) tj^tjbticd) + ^(abIcd)) dV (14) 

(15) 



Similarly, as the equilibrium configuration is symmetric with re- 
spect to the equatorial plane, the integrations are performed over 
symmetric intervals in the vertical direction and the coupling 
coefficients vanish when the integrands are odd functions of z. 
This happens when an odd number of modes is involved, whose 
Lagrangian eigenfunctions are antisymmetric with respect to the 
reflection z «-» -z, 



\A\ 



Y\e Ak =-1 => K A = 0. 



(22) 



1 



In the above equation, ex s is the parity of the A^.-mode under re- 
flection about the equatorial plane, and = ±1 corresponds to 
even and odd modes, respectively. The condition d22l is referred 
to as the vertical selection rule. 



3.2. Nodal modes 

We consider modes whose Lagrangian eigenfunctions have at 
least one (but still less than 1 //?) node in both the r and z di- 
rections; such modes are referred to as the nodal modes. The 
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Table 1. Selected coupling coefficients of the epicyclic modes. 
The second column provides formulae for the coupling coeffi- 
cient in a general gravitational field. The third column provides 
the corresponding numerical values of the dimensionless cou- 
pling coefficients for the case of the pseudo-Newtonian potential 
(see section 2]i. The radius of the torus is set to the location of 
the 3:2 epicyclic resonance, ro = 



/(•(general formula) 



* (PN) 



^rvv — "f vv 



4^ [^2 <Pl2 + <P2l] 



-V2 
18.7539 

-11.6211 

-112.754 



Lagrangian displacements and their gradients satisfy the scal- 
ings 

£a~Wa, V£a ~ ft A /(J3ro). (23) 
The coupling coefficients can be estimated by 

xf ~ Po V [M/(J3r )] w , xf ~ poV<S> [M/r ] w , (24) 

where the torus volume and the central pressure are 

V~p 2 rl po ~ p 2 r 2 Q Q. 2 p (25) 

Therefore, the gravitational coupling coefficients are smaller 
than the pressure ones by a factor proportional to 



(26) 



Nonlinear dynamics of the nodal modes is therefore dominated 
by pressure coupling (we note that \A\ > 3). 

3.3. Epicyclic modes 

The Lagrangian eigenfunctions of the epicyclic modes are al- 
most uniform on the torus cross-section, They cannot be cou- 
pled with each other by pressure forces in the infinitely slender 
torus, since they are determined by gradient s in the Lagrangia n 
displacements. However, as pointed out by Blaes et al. (2007), 
they may be coupled by a non-slender corrections in somewhat 
thicker tori. 

To examine this possibility we approximate the eigenfunc- 
tions of the epicyclic modes to be 



^+/3^+0(J3 2 ), 



(27) 



where ^ are given by (0 and (0, and both coefficients are 
roughly similar, \§P V \ ~ J[ uw . We calculate the leading orders for 
both the pressure and the gravitational coupling coefficients. The 
contribution of the zeroth-order eigenfunctions to the gradient 
in the Lagrangian displacement is highly reduced because the 
eigenfunctions are strictly uniform over the torus cross-section 
to this order of approximation. The only nontrivial corrections to 
the gradient are due to a small azimuthal curvature of the torus; 
they are however of the same order in (3 as the contributions of 
the non-slender corrections. The pressure coupling coefficients 



are therefore even smaller than the gravitational coupling coeffi- 
cients by a factor of 



(28) 



and the total coupling coefficients are scaled according to the 
gravitational contributions, 



(29) 



A -"•"0'0' 

We conclude that in slender tori the nonlinear coupling between 
the epicyclic modes is governed by gravity. 

3.4. Axisymmetric epicyclic modes 



As demonstrated in Sect. 12.11 each epicyclic mode of the torus 
is accompanied by its "complex-conjugated" analog, whose 
frequency is opposite and whose eigenfunction is complex- 
conjugated. It is therefore useful to introduce the notation 



-0>A, 



£a - hA 



(30) 



and denote complex-conjugated pairs by (A, A). 

The second column of Table Q] contains few examples of 
the coupling coefficients between the axisymmetric epicyclic 
modes obtained by substitution of relations © and ([3j into equa- 
tions ( fT3l . ( TT3T > and (TTTb - The integrands were approximated by 
leading-order terms in their /3-expansions. For a simpler nota- 
tion, we introduce (py = (<9j.<9j<E>)o. 

4. Epicyclic resonances 

We examine nonlinear interaction between two epicyclic modes 
of the torus. We assume that the resulting oscillations can be 
described by the Lagrangian displacement of the form 

£(?, x) = c r (f) + a (t) + c v (t) f v (x) + cv(?) &(*)■ (3D 

Any influence of the other modes on the dynamics of these two 
is ignored (we return to this issue in Sect. 15.3b . The reality of 
the Lagrangian displacement requires that c f = c* and Cv = c*. 
The coefficients c r (f) and c v (f) are solutions of the two coupled- 
oscillator equations ©. We assume that ca ~ e <k 1 (e is a 
dimensionless parameter, describing the smallness of the ampli- 
tudes), then the oscillators are weakly coupled and the nonlin- 
earities on the right-hand sides of the equations can be treated as 
perturbations. 

Apart from the main oscillations, whose frequencies are 
close to the eigenfrequencies o) t = m r Qo + to r and ai v = m Y Q.o + 
oj z , the solutions contain harmonics that appear because of the 
nonlinearities in equations ([6j. When the two eigenfrequencies 
are almost commensurable, interactions between the harmonics 
and main oscillation lead to resonances. The order of resonance 
is given by the order of nonlinearities required to excite the res- 
onance. For a general system with two degrees of freedom, the 
order of resonance p:q (p and q are relative prime numbers) is 
given by 



p + q-l. 



(32) 



The only exception to this simple rule is the 1 : 1 resonance, the 
order of which is = 3. 

Generally, the strength of the resonance decreases with in- 
creasing order. The resonances of higher order are more diffi- 
cult to tune because the resonance range scales as e"«o» (a> is 
an eigenfrequency of the system). Moreover, the amplitudes and 
phases of resonant oscillations are modulated on the timescale 
proportional to e"'' «<y _1 . 
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Table 2. Possible resonances up to the fourth order. Due to the 
equatorial -plane reflection symmetry of the equilibrium torus, 
many resonances are absent between epicyclic modes. 

Order General system Epicyclic modes 

resonance resonance condition 



2nd 


1:2 


1:2 


m x = 2m v 




2:1 






3rd 


1:3 








1:1 


1:1 


ra r = ra v 




3:1 






4th 


1:4 


1:4 


m r = 4m v 




2:3 








3:2 


3:2 


3ra r = 2m v 




4:1 







4.1. Selection rules for epicyclic resonances 

In systems with some intrinsic symmetries, the presence of har- 
monics and occurrence of resonances depends on the symmetry 
properties of oscillation modes. Some resonances do not occur 
even though the eigenfunctions satisfy the corresponding reso- 
nance conditions because the nonlinearities involved in the pro- 
duction of the harmonics vanish. 

We explore possible resona nces up to the fourth or der using 
the method of multiple scales (Nayfeh & Mook 1979I) . The re- 
sult is shown in the second column of Table[2] We assume that all 
coupling coefficients are nonzero, which corresponds to a gen- 
eral system with no intrinsic symmetry. Next, we apply the ver- 
tical selection rule ( l22b . taking into account parities of the radial 
and vertical epicyclic modes, e l v = + 1 . Possible resonances are 
listed in the third column of Table [2] The fourth column shows 
additional conditions for azimuthal wavenumbers obtained using 
the azimuthal selection rule ( 1211 ). 

Apparently, both selection rules reduce significantly the 
number of possible resonances - more of them are possible only 
for the axisymmetric modes (m r = m w = 0). In the strong grav- 
itational field of both rotating and non-rotating black holes, the 
vertical epicyclic frequency is always greater than the radial one. 
Therefore the first three resonances listed in the third column of 
Table [2] does not occur for axisymmetric epicyclic modes and 
the strongest epicyclic re s onanc e is therefore 3:2. We note that 
iKluzniak & Abramowic^ d2002l) anticipated this result by using 
the analogy of t he pa rametric resonance in t he Mathieu equation. 
Rebusco (2004]) and lHorak & Karasl ((2006) achieved similar re- 
sults in their discussion of internal resonances in the test-particle 
epicyclic motion. 

As follows from the azimuthal selection rule, the neces- 
sary condition for p:q resonance between two non-axisymmetric 
epicyclic modes is pm r - qm v = (see Table [2J. Therefore, 
the resonance condition is the same as for the axisymmetric 
modes, (o z /o> r = p/q (and the resonance occurs therefore at 
the same radius), but the frequencies of the modes are different, 
u> v — m v D. ± o> : and o> x — m x Q. + co r . 



4.2. The strong-gravity 3:2 epicyclic resonance 

We explore the strength and resonance range of the 3:2 
epicyclic resonance between the radial and vertical axisymmet- 
ric epicyclic modes in a slender torus. We study parametric exci- 



tation of vertical epicyclic oscillations by radial oscillations. For 
simplicity, we ignore any feedback of the vertical oscillations to 
the radial. This is a reasonable approximation when the ampli- 
tude of the radial mode is far greater then the amplitude of the 
vertical mode. 

The effects of strong gravity on the central object (such as a 
non-rotating bl ack hole or a compact n eutron star) are included 
by using the iPaczvriski & Wiital dl980 ) pseudo-Newtonian po- 
tential, 



GM 
R-R, 



R 



V'" 2 + z 2 , 



(33) 



where G is the gravitational constant, M is the mass of the com- 
pact object and R s = 2GM/c 2 is the Schwarzschild radius. The 
radial epicyclic frequency is always smaller than the vertical one. 
The latter is equal to the Keplerian orbital frequency because the 
gravitational field of the compact object is spherically symmet- 
ric. The 3:2 epicyclic resonance occurs when the torus radius 
is r = r 3:2 = 9.2GM/c 2 (the 3:2 radius is at 10.8 GM/c 2 in 
Schwarzschild spacetime). 

We first renormalize the coefficients ca in Eq. © and intro- 
duce the dimensionless coupling coefficients ka, 



ca 



u) A b A 



y At 



-C A , K A 



0"() 



< A < [Mrffi] 



n 

k=l 



The governing equations then take the form 
— - + \u> A c A = uo A T A {ci), 



(34) 



(35) 



where the formula for the dimensionless nonlinear functions 
Ta(ci) is analogous to Eq. ([8]) (all quantities are indicated by 
a bar). Numerical values of several dimensionless coupling co- 
efficients evaluated at ro = r^a are listed in the third column of 
Table □ 

The solution of equation d35l l can be found using the method 
of multiple scales. The radial oscillations can be approximated 
by c r = A\ exp(-iw,.f). Without any loss of generality, we as- 
sume that the initial phase of the oscillations is such that A t e R. 
Up to the leading order, the vertical oscillations are given by 
c v = A v (f) exp(-iw.?X where A v (f) € C is a slowly changing 
amplitude. The slow time evolution is given by the amplitude 
equation, the form of which is 



dA v 



6t 



(^r|A r 



■A V \A V \ 2 )A V +«A;I I 3 e ifa ' = 0, 



(36) 



where 8o> = 2u> z - 3u r describes the detuning of the eigenfre- 
quencies from the rational ratio. It depends on the location of 
the torus center ro (6o> = when ro = r^i). The dimensionless 
constants a, A r and A v are given by coupling coefficients. The 
corresponding formulae are 



1 

— Cxi- 

2 1 



9^rvv(^rrvv "I" ^itvv) 8/£*n 



3 



(37) 

(38) 
(39) 



Numerical values of these constants at the resonance radius are 
a = 1101, A T = 78.14 and A v = 1.333. 
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Fig. 1. Region of the 3:2 parametric-like resonance between the 
epicyclic modes in pseudo-Newtonian slender tori. The reso- 
nance region is shown in the plane of the torus position (horizon- 
tal) versus the radial epicyclic mode amplitude (vertical). The 
domains where the vertical epicyclic oscillations are stable and 
unstable are separated by transition curves. The contours denote 
constant growth-rate of vertical oscillations. The corresponding 
values of y/Qo are multiples of 0.25 x 1CT 2 (zero corresponds to 
the transition curves). The inset of the figure shows the position 
and radial amplitude of the torus ("STA"-point) in the simula- 
tions of Sram kova etail d2007l) . 



When the amplitude of the vertical oscillations is small, the 
second term in brackets in Eq. (l36l l can be neglected with respect 
to the first one and we obtain a linear equation. We attempt to 
find its solution in the form 



A v = [«!(?) +ia 2 (?)]e 



-iSa)t/2 



(40) 



with two unknown real functions a\(i) and a^if). Substituting 
these expressions into equation (l36l >. we obtain 



doi 
~df 

da 2 

~d7 



+ [(5« + cj z (A T Af + aAfj] a 2 = 0, 
+ [dcj + u z {ArAj - aAfj] a\ = 0. 



(41) 
(42) 



These linear equations have constant coefficients, therefore, we 
may assume that a\ i depends on / as a\2 - 5i,2 exp(-yf)- When 
this form is inserted into equations ( l4TT i and (l42l we obtain a 



couple of linear homogenous algebraic equations for a 1,2 
are solvable when 



that 



r = 



a 2 o? z A\ - ( ^6a> + lj v A x A 2 



1/2 



(43) 



The amplitude of the vertical oscillation grows exponentially 
when the frequency detuning is in the range u)\ < 6 to < ci>2 
with 



(44) 



The region of resonance is shown in Fig.Q] The transition curves 
separating the stable and unstable regions are given by Eq. d44l . 
For a given amplitude of the radial oscillations, the size of the 
resonance range Aw = \<x>2 - wil is proportional to A]. Instead 
of the frequency detuning 5u, we use the corresponding radial 
coordinate ro of the torus center in the figure (ro is calculated 
from the well-known dependences of the epicyclic frequencies 
on the radius). 

The growth-rate of the unstable vertical mode is given by 
Eq. ( l43l and indicated by contours in Fig. Q] The maximal 
growth-rate that can be obtained for a given amplitude of ra- 
dial oscillations is proportional to A;? and occurs when 5co = 
—2to z A T A 2 . It is given by y max = au) z A]. In the figure, the am- 
plitude of the radial epicyclic oscillations is shown in the units 
of GM/c 2 . The value of A r is connected to the value A r of the 
dimensionless amplitude by 



2£ 



Mr 2 Q 2 _ V2r 

A- — z A r , 



(45) 



where ^[ is the radial component of the radial epicyclic-mode 
eigenfunction [Eq. @]. 



5. Discussion 

5. 1 . Strength of the resonance 

We have found that the 3:2 epicyclic resonance is very sensitive 
to the precise tuning of the eigenfrequencies of the torus. For 
small amplitudes of the radial oscillations, Eq. (l44l > implies that 
the range of the detuning parameter, for which the resonance 
may operate, is very limited. Moreover, the growth-rate of the 
vertical epicyclic oscillations in the resonance is quite small. It 
appears that these facts make the epicyclic resonance difficult to 
observe in both, numerical simulations as well as real astrophys- 
ical objects. 

In principle, our results agree with the numerical simula- 
tions of Sram kova et alj d2007l) . who claim that the epicyclic 
modes are not resonantly coupled. In their simulations the ra- 
dial extent of the torus is smaller than its radius by the factor 
~ 0.02, while the velocity amplitude of the radial perturbations 
are smaller than the central sound speed by the factor ~ 0.3, 
This situation corresponds to the amplitude of the radial oscil- 
lations A r ~ 10~ 2 GM/c 2 . In the simulations, the torus center 
is at ro = 9.2GM/c 2 (models A3 and A4). The corresponding 
point [ro, A r ] is outside the resonance tongue as shown in Fig.Q] 
For the same amplitude of radial oscillations our theory predicts 
maximal growth-rate y max ~ 10~ 7 Oo- Besides the high precision 
of the calculations, an observation of the resonance in numeri- 
cal simulations requires considerably long simulation times and 
a precisely tuned radius of the torus. 
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5.2. Comparison of the characteristic time-scales 

In contrast to the inviscid flow considered in this paper, real ac- 
cretion tori are made of viscous fluids. Low viscosity causes 
slow secular evolution of the torus by changing its angular- 
momentum distribution on a viscous timescale. These changes 
affect the torus radius ry. Similarly, a viscous heating causes sec- 
ular evolution of the torus volume on the thermal timescale. Both 
processes influence the eigenfrequencies of the torus and may 
therefore disturb the precise tuning, which is required by the res- 
onance. This happens when the resonant modulation timescale 
l/y exceeds the characteristic time that the torus spends in the 
resonance region. We adopt common formulae for a thin disk 
(iFrank et alJll992l) 

f visc = a; 1 (H/R)- 2 forb ~ a; l /3~ 2 t olb , (46) 

fth = o^'foi-b, (47) 

where a t parameterize the accretion disk turbulence 
dShakura & Svunvaevl 1 1973b and f 01 b ~ C2 ' is the orbital 
timescale. The typical drifts of eigenfrequencies on these 
timescales due to the changing radius or the volume of the 
torus are (Aw) v ; sc ~ Qo and (Aco)^ ~ /3 2 Qo. Therefore, 
typical times that the torus spends in the resonance re- 
gion undergoing viscous or thermal changes are similar, 
(Ar) visc ~ (Ar) th ~ or" 1 /r 2 a(A r /ro) 3 f or b- If we compare this 
result with the resonant timescale, we obtain an upper limit on 
a t for the resonance to operate, 

(3 2 a t <a 2 (^j. (48) 

5.3. Effects of other oscillation modes 

We considered nonlinear interactions between two epicyclic 
modes and ignored the influences of all other modes [see 
Eq. d3"Ttl. If they are present with small amplitudes in the os- 
cillations, the Eq. (t36l > is modified by the presence of additional 
terms /l/|A/| 2 in brackets. The shape of the transition curves and 
the growth-rates of the vertical epicyclic modes are then slightly 
changed, however the size of the resonance range and the general 
discussion presented in this paper remain the same. 

More importantly, the epicyclic resonance may be sup- 
pressed by the parametric instability that may quickly advect 
the energy from t he epicyclic modes to some low-frequency 
modes of the torus (lDziembowskilll982l:IWu & GoldreicM200U 
Arra s et ail 120031: fNowakowski 2005). Being a resonant inter- 
action among three modes, the characteristic timescale of the 
parametric instability is far shorter than that of the 3:2 epicyclic 
resonance. An important condition for this process is the ex- 
istence of a pair of low-frequency modes that form a resonant 
triple with the epicyclic modes. This happens when the frequen- 
cies and azimuthal wavenumbers of the low-frequency modes 
satisfy the conditions Lu> rz « (±>\ + u>2 and m\ + ni2 = 0. The 
'parent' mode (here radial or vertical axisymmetric epicyclic 
mode) of the highest frequency u> nz is a source of energy for 
the two 'daughter' modes with frequencies 0*1,2. For each pair 
of daughter modes, there exists a lower limit to the parent-mode 
amplitude above which the parametric instability begins to oper- 
ate. This limit depends on the damping rates j\, j2 of daughter 
modes and on a couplin g coefficient k of the three-mode interac- 
tion as A m i n oc 7172/kl 2 (Dziembowski 1982). 

To decide whether the parametric instability advects energy 
from the epicyclic modes to some other low-frequency (per- 
haps unobserved) modes, it is necessary to explore in details the 



eigenfrequencies, eigenfunctions and damping rates of the torus 
modes . Su ch analyses has be en carried out by I Wu & Goldreichl 
d200ll) and lArraset a l. (2003) in the context of stellar pulsations. 
A similar study is beyond the scope of this paper because the 
damping processes in tori are still only purely understood. We 
therefore summarize the necessary properties of potential daugh- 
ter modes, based on known properties of the eigenfrequencies 
and eigenfunctions of the slender torus. 

We first examine whether possible pairs of daughter modes 
exist among the lowest orde r modes of the c onstant angular- 
momentum tori derived by iBlaes et ail d2006l) (see their Table 
4). Both axisymmetric and non-axisymmetric modes with m 
I and 2 are considered (the modes with higher azimuthal 
wavenumber are excluded because their frequencies are higher 
than (L>r, z )- We evaluate the eigenfrequencies in the pseudo- 
Newtonian gravitational field at r — ry^ and select only those 
pairs whose mutual interactions with the epicyclic modes are not 
forbidden by the selection rules. For each such combination, we 
evaluate the frequency detunings Au> nz = cj v - <±>\ - ci>2- We 
do not find any pair that provides Au> rtZ within 10% of the or- 
bital velocity Qo. The situation is probably similar for steeper 
angular-momentum distributions. We expect that the resonance 
condition will be satisfied only for some particular values of k. 
Therefore, we may conclude that the resonant coupling is inef- 
ficient and the parametric instability does not drive the energy 
from the epicyclic to the lowest-order modes of the torus. 

The high-order modes, whose eigenfunctions have many 
nodes in both the r and z directions, can be treated using the 
WKBJ approximation. These modes are damped more strongly 
and therefore they advect energy from the epicyclic modes more 
effectively. Their eigenfunctions can be approximated by 

£ = | e i(»«A+^ # = j krdr + j fezdZ) (49) 

where k r ,k z » rr and t- are slowly varying functions of r and 
z and $ is a rapidly changing phase. The wave-vector (k r , k z ) is 
connected to the eigenfrequency of the mode by the dispersion 
relation that takes the form (i) <x 2 = c 2 (k 2 + k 2 ) for the acoustic 
and surface-gravity modes (c s is local speed of soun d) and (ii) 
cr 2 = [k 2 /(k 2 + k 2 )]k 2 for the internal-inertial modes (Blae s et ail 
I2006T) . 

We suppose that the two high-order daughter modes, 1 and 
2, are governed by the same dispersion relation. The integrands 
of the three-mode coupling coefficients /e r i2 and K v i2 contain an 
exponential exp[i(#i + §2)] (the epicyclic modes are almost uni- 
form, therefore their contributions are negligible). The integrand 
contributes significantly only close to the point where the phase 
of the exponential becomes stationary in both r and z directions 
and therefore k\ r + k2 r ~ k\ z + k2 Z ~ 0. Using the dispersion 
relation, we conclude that this happens when cr\ m ±0-2- On 
the other hand, the two modes form a resonant triple with an 
epicyclic mode when oi rA su cr\ + (r%. Hence, the corotation fre- 
quencies of the potential daughter modes are cr\ « <x 2 as u> nz /2, 
which corresponds to the eigenfrequencies cjj^ ~ <w r , z /2 + mQ(>. 
The parametric instability operates only when the mode offer- 
ing the energy is the one with the highest frequency, i.e. when 
wi,2 < <L>r,z- The only possible value of the azimuthal wavenum- 
ber is then m = 0. 

The frequencies of the acoustic and surface-gravity modes 
that are governed by the dispersion relation (i), increase with in- 
creasing k r and k z . Therefore these modes will not drain energy 
from the epicyclic modes. Consequently, the three-mode para- 
metric instability is not dangerous for the epicyclic resonance in 
constant angular-momentum tori. 
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The frequencies of the inertial oscillations governed by rela- 
tion (ii) are, however, always between and k. As the angular- 
momentum distribution approaches the Keplerian one, k ap- 
proaches the radial epicyclic frequency. Therefore, for suffi- 
ciently steep distribution of angular momentum, there will be 
high-order axisymmetric inertial modes whose eigenfrequencies 
are sufficiently close to a» rjZ /2. These modes may be important 
for the parametric instability. Further careful analysis is required 
in order to determine whether the nonlinear interactions with 
these modes is sufficient to suppress the epicyclic resonance. 

6. Conclusions 

We have applied a general theory of nonlinear pulsation in rotat- 
ing stars to a problem of nonlinear oscillations of thick accretion 
disks. In this note, we have calculated the strength of the cou- 
pling between epicyclic modes of slender torus. We have taken a 
closer look to a parametric excitation of vertical epicyclic motion 
of the torus due to radial epicyclic oscillations in the epicyclic 
resonance. Our main findings can be summarized as follows. 

(1) While the nodal modes of the torus are coupled by pres- 
sure gradients, the coupling between the epicyclic modes is gov- 
erned by gravity (see Sect. [3). 

(2) The strongest epicyclic resonance between the axisym- 
metric modes occurs when the ratio of the vertical and hori- 
zontal epicyclic frequencies are close to the 3:2 ratio. In ad- 
dition, the 3:2 resonance may appear also between two non- 
axisymmetric modes, whose azimuthal wavenumbers satisfy 
3m r = 2m v . However, already the first possible resonance of this 
kind (m r = 2 and m v = 3) provides frequencies that are ap- 
parently too high to be identified with those observed in QPOs 
(Sect.©. 

(3) The azimuthal selection rule forbids any resonance be- 
tween one axisymmetric and one non-axisymmetric epicyclic 
mod e. No t e that such resonance has been recently suggested 
by iBursal d2005l) to explain discrepancy between spin esti- 
mates based on spectral fitting and resonance models. We note 
that the conclusions (3) and (4) are consequences of the ax- 
ial and equatorial-plane symmetry of the equilibrium torus only 
(Sect.©. 

(4) We have examined the stability of the vertical epicyclic 
motion in the 3:2 epicyclic resonance with the radial epicyclic 
oscillations. For a given amplitude of the radial oscillations we 
have found a range of torus radii for which the vertical oscilla- 
tions are unstable. We have noted that the epicyclic resonance 
become significant only for high amplitudes of the radial oscil- 
lations (see Sect. 14.21 ). 

Acknowledgements. I appreciate fruitful discussions with Paola Rebusco, Eva 
Sramkova, Marek Abramowicz, Wlodek Kluzniak, Omer Blaes, Gabo Torok 
and colleagues from the Astronomical institute in Prague. I am also grateful 
to an anonymous referee for many valuable comments that largely improved the 
manuscript. Finally, I acknowledge the hospitality of MPI Garching and finan- 
cial support of the GACR grant 205/06/P415 and of the Center for Theoretical 
Astrophysics (LC 06014). 



Blaes O.M., Arras P., Fragile PC, 2006, MNRAS, 369, 1235 
Blaes O.M., Sramkova E., Abramowicz M.A., Kluzniak W., Torkelsson U., 2007, 
ApJ, 665, 642 

Bursa M., 2005, In: Hledfk S., Stuchlfk Z. (eds.) RAGtime 6/7: Workshops on 

black holes and neutron stars, 39^5 
De Villiers J.P., Hawley J.F., Krolik J.H., 2003, ApJ, 599, 1238 
Dyson J., Schutz B.F., 1979, Royal Society of London Proceedings Series A, 

368,389 

Dziembowski W., 1982, Acta Astronomica, 32, 147 

Fogelstrom S., Levin L., Slapak R., 2008, ArXiv:0802.1594 

Frank J., King A., Raine D., 1992, Accretion Power in Astrophysics, Cambridge 

University Press 
Goldreich P., Goodman J., Narayan R., 1986, MNRAS, 221, 339 
Goodman J., Narayan R., Goldreich P., 1987, MNRAS, 225, 695 
Horak J., Karas V., 2006, A&A, 451, 377 
Kato S., 2003, PASJ, 55, 801 
Kato S., 2004, PASJ, 56, 905 
Kato S., 2008, PASJ, 60, 111 

Kluzniak W., Abramowicz M.A., 2002, astro-ph/0203314 

Kumar P., Goldreich P., 1989, ApJ, 342, 558 

Lee W.H., 2005, Astronomische Nachrichten, 326, 838 

Lee W.H., Abramowicz M.A., Kluzniak W., 2004, ApJ, 603, L93 

Machida M., Nakamura K.E., Matsumoto R., 2006, PASJ, 58, 193 

McClintock J.E., Remillard R.A., 2006, In: Lewin W., van der Klis M. (eds.) 

Compact stellar X-ray sources, 157-213, Cambridge University Press 
Narayan R., Goldreich P., Goodman J., 1987, MNRAS, 228, 1 
Nayfeh A.H., Mook D.T., 1979, Nonlinear Ocillations, Wiley-Interscience, New 

York 

Nowakowski R.M., 2005, Acta Astronomica, 55, 1 
Paczynski B., Wiita P.J., 1980, A&A, 88, 23 
Papaloizou J.C.B., Pringle J.E., 1984, MNRAS, 208, 721 
Papaloizou J.C.B., Pringle J.E., 1985, MNRAS, 213, 799 
Rebusco P., 2004, PASJ, 56, 553 

Rezzolla L., Yoshida S., Maccarone T.J., Zanotti O., 2003, MNRAS, 344, L37 
Schenk A.K., Arras P., Flanagan E.E., Teukolsky S.A., Wasserman I., 2002, 

Phys. Rev. D, 65, 024001 
Schutz B.F., 1980a, MNRAS, 190, 7 
Schutz B.F., 1980b, MNRAS, 190, 21 
Shakura N.I., Syunyaev R.A., 1973, A&A, 24, 337 
Sramkova E., Torkelsson U, Abramowicz M.A., 2007, A&A, 467, 641 
Van Hoolst T., 1994, A&A, 286, 879 
Van Hoolst T., Smeyers P., 1993, A&A, 279, 417 

Vio R., Rebusco P., Andreani P., Madsen H., Overgaard R.V., 2006, A&A, 452, 
383 

Wu Y, Goldreich P., 2001, ApJ, 546, 469 



References 

Abramowicz M.A., Kluzniak W., 2001, A&A, 374, L19 

Abramowicz M.A., Karas V, Kluzniak W., Lee W.H., Rebusco P., 2003, PASJ, 
55, 467 

Abramowicz M.A., Blaes O.M., Horak J., Kluzniak W., Rebusco P., 2006, 

Classical and Quantum Gravity, 23, 1689 
Arras P., Flanagan E.E., Morsink S.M., et al., 2003, ApJ, 591, 1 129 
Balbus S.A., Hawley J.F., 2002, ApJ, 573, 749 
Blaes O.M., 1985, MNRAS, 216, 553 



